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1 .  Introduction 


As  a  result  of  research  conducted  over  the  past  20  years,  the  capability 
now  exists  to  unambiguously  identify  almost  all  seismic  events  having 
magnitudes  characteristic  of  well-coupled,  underground  nuclear  explosions 
with  yields  greater  than  a  few  kilotons  (i.e.  mb  >  4,  OTA  (1988)).  However, 
in  the  context  of  monitoring  any  eventual  Comprehensive  Test  Ban  Treaty 
(CTBT),  consideration  will  have  to  be  given  to  the  ci^ability  for  identifying 
the  much  smaller  signals  which  might  be  expected  to  result  from  possible 
evasive  testing  involving  cavity  decoupling.  That  is,  since  the  U.S.  nuclear 
cavity  decoupling  experiment  STERLING  has  established  that  it  is  possible 
to  reduce  the  aiiq)litude  of  the  radiated  seismic  signal  by  a  least  a  factor  of 
70  using  this  testing  procedure,  it  follows  that  cominehensive  monitoring  of 
underground  nuclear  tests  in  the  1  to  10  kt  range  will  necessarily  involve 
identification  analyses  of  small  seismic  events  with  magnitudes  in  the  range  2 
<  mb  <  3.  Moreover,  at  such  low  magnitudes,  naturally  occurring 
earthquakes  are  supplemented  by  a  large  number  of  chemical  explosions 
(CE)  of  similar  size  which  are  associated  with  the  variety  of  quarrying, 
mining  and  construction  projects  which  are  continuously  being  carried  out  in 
most  developed  areas  of  the  world.  This  constitutes  a  potentially  serious 
monitoring  problem  in  that  there  are  no  consistently  reliable  procedures 
available  at  the  present  time  for  distinguishing  between  seismic  signals 
produced  by  cavity  decoupled  nuclear  and  chemical  explosions.  Thus,  even 
assuming  the  existence  of  the  very  extensive  monitoring  networks  needed  to 
detect  such  small  events,  questions  remain  concerning  our  capability  to 
correctly  identify  the  signals  from  small  decoupled  explosions  among  the 
numerous  signals  expected  from  earthquakes  and  CE  events  of  comparable 
size.  The  objective  of  the  research  program  described  in  this  report  is  to 
improve  the  capability  to  distinguish  between  small  cavity  decoupled  nuclear 
explosions  and  CE  events  through  analyses  of  simulated  and  observed 
seismic  data  representative  of  the  two  source  types. 
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This  report  presents  a  summary  of  the  preliminary  investigations 
which  have  been  conducted  during  the  first  year  of  this  project.  In  Section 
2,  near-regional  seismic  data  recorded  from  the  U.S.  nuclear  cavity 
decoupling  test  STERLING  and  the  nearby,  tamped  CE  test,  denoted 
STERLING  HE,  are  analyzed  and  compared  in  an  attempt  to  identify  any 
diagnostic  differences  between  these  two  source  types.  This  is  followed  in 
Section  3  by  an  analysis  in  which  seismic  data  recorded  at  NORSAR  from 
Soviet  nuclear  tests  conducted  in  salt  cavities  at  the  Azgir  site  north  of  the 
Caspian  Sea  are  compared  with  corresponding  data  recorded  at  the  same 
stations  from  events  of  about  the  same  size  which  appear  to  be  CE  events 
conducted  in  the  vicinity  of  the  Azgir  site.  In  Section  4,  broadband  seismic 
data  recorded  at  near-regional  Soviet  stations  from  an  Azgir  cavity 
decoupling  test  are  theoretically  scaled  and  compared  with  data  recorded  at 
similar  distances  at  the  nearby  IRIS  sudon  KIV  from  well-documented  CE 
events.  The  report  concludes  with  Section  5  which  contains  a  summary  and 
statement  of  preliminary  conclusions  regarding  the  identification  of  seismic 
signals  recorded  from  cavity  decoupled  nuclear  and  chemical  explosions. 
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2.  Comparison  of  Broadband  Near-Regional  Seismic  Signals 
Recorded  from  the  STERLING  and  STERLING  HE 
Explosions 


The  SALMON,  STERLING  and  STERLING  HE  explosions  were 
conducted  in  the  Tatum  salt  dome  near  Hattiesburg,  Mississippi  between 
1964  and  1966  with  the  objective  of  evaluating  the  effectiveness  of  cavity 
decoupling  in  reducing  the  amplitude  of  the  radiated  seismic  signal  from 
underground  nuclear  explosions.  The  first  of  these  tests  was  SALMON, 
which  was  a  fully  coupled  5.3  kt  nuclear  explosion  conducted  at  a  depth  of 
828  m  on  October  22,  1964.  This  was  followed  on  November  17,  1966  by 
STERLING  HE,  a  2.7  ton  tamped,  high  explosive  calibration  shot  which  was 
detonated  about  350  m  southwest  of  SALMON  at  a  depth  of  831  m.  The 
decoupled  nuclear  test  STERLING  had  a  yield  of  0.38  kt  and  was  detonated 
on  December  3,  1966  in  the  17  m  radius  semispherical  cavity  produced  by 
the  SALMON  explosion  (Ferret  1968a,b;  Springer  et  al.  1968).  Given  that 
the  low  frequency  decoupling  factor  for  STERLING  has  been  estimated  to 
have  been  q)ixoximately  70  ±  20  (Springer  et  al.,  1968),  it  follows  that  the 
an[q>litudes  of  the  low  frequency  seismic  signals  from  STERLING  were  of 
the  same  order  of  magnitude  as  those  to  be  expected  from  a  5.4  ton  tamped 
nuclear  explosion.  It  follows  that  since  HE  is  generally  thought  to  produce 
low  frequency  seismic  signals  which  are  about  twice  as  large  as  those 
produced  by  a  nuclear  explosion  of  con^arable  energy  release,  the  observed 
seismic  data  from  STERLING  should  be  directly  comparable  to  that 
observed  from  the  2.7  ton  STERLING  HE  explosion.  This  hypothesis  is 
supported  by  the  fact  that  the  observed  seismic  amplitude  levels  at  fixed 
recording  stations  were  generally  quite  similar  for  these  two  events.  Thus, 
these  data  provide  a  nearly  ideal  opportunity  to  compare  seismic  data 
recorded  from  tamped  CE  and  nuclear  cavity  decoupled  explosions  which 
are  nearly  co-located  and  have  very  similar  seismic  source  strengths. 
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Although  data  from  STERLING  and  STERLING  HE  were  recorded  at 
a  number  of  stations  located  within  about  70  km  of  their  common  epicenter, 
the  miyority  of  these  data  are  of  questionable  reliability  due  to  signal-to- 
noise  limitations  and  other  problems.  Therefore,  in  the  present  analysis  we 
will  focus  on  the  stations  denoted  as  lOS  and  20S  which  were  located  at 
ranges  of  16  and  32  km  south  of  the  source  locations,  respectively  (Murphy, 
1969).  Figures  1  and  2  show  comparisons  of  the  vertical  component, 
broadband  (1-35  Hz)  recordings  from  STERLING  and  STERLING  HE  at 
stations  lOS  and  20S.  It  should  be  noted  that  although  only  these  vertical 
component  data  will  be  considered  in  the  following  discussion,  it  has  been 
found  that  the  other  components  of  motion  display  similar  characteristics 
and,  therefore,  the  conclusions  should  be  considered  to  be  representative  of 
them  as  well.  Figure  1  shows  8  second  segments  of  the  recorded  time 
histories  at  station  lOS  from  STERLING  (top)  and  STERLING  HE 
(bottom).  These  records  are  relatively  complex  for  both  sources  at  this 
near-regional  distance,  presumably  as  a  result  of  the  propagation  path 
heterogeneities  introduced  by  the  source  region  salt  dome  structure  and  the 
surrounding  sedimentary  layered  sequence  which  it  pierces.  Many  consistent 
phases  are  apparent  in  the  records  between  the  two  events.  The  initial  P- 
wave  windows  appear  to  match  quite  closely  fcr  the  first  second  or  so. 
Beginning  roughly  two  seconds  after  the  first  arrival  is  a  second  phase 
consisting  of  a  relatively  low-frequency  carrier  with  high-frequency  signals 
superinq>osed.  This  latter  phase  produces  the  largest  amplitudes  on  the 
vertical  component  records,  and  its  interpretation  is  not  clear.  Springer  et 
al.  (1968)  identified  this  signal  as  an  S  phase,  but  Gupta  et  al  (1987) 
determined  on  the  basis  of  synthetic  seismograms  derived  for  a 
representative  structure  that  this  phase  arrived  too  early  for  S  and  was  more 
likely  reverberating  P-waves  in  the  sediments.  Accepting  this  latter 
interpretation  the  S-wave  window  begins  about  one  second  later;  we  follow 
this  interpretation  throughout  the  subsequent  discussion.  Comparing  the 
signals  from  the  two  source  types,  it  can  be  seen  that  the  most  prominent 
difference  is  in  the  ratio  of  the  S  to  initial  P  amplitudes,  which  appears  to  be 
on  the  order  of  a  factor  of  three  larger  for  STERLING  HE  than  for 
STERLING. 
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Similarly,  Figure  2  shows  16  second  segments  of  the  ground  motion 
time  histories  recorded  at  station  20S  from  STERLING  (top)  and 
STERLING  HE  (bottom).  The  records  are  again  fairly  complex  including 
multiple  arrivals  in  the  P-wave  window  which  are  iy>parently  associated  witii 
the  propagation  path  since  we  observe  general  consistency  between  the  two 
sources.  In  this  case  the  largest  observed  amplitudes  begin  arriving  about 
five  seconds  after  the  initial  P.  This  phase  is  generally  lower  in  frequency 
and  is  thought  to  be  S;  in  this  case  Gupta  et  al.  (1987)  agree  with  the 
Springer  et  al.  (1968)  interpretation.  Just  as  for  P,  the  S  arrival  has  a  long 
duration  with  multiple  arrivals  i^parently  associated  witii  propagation  in  the 
complex  sedimentary  structure.  Once  again,  the  most  prominent  difference 
between  these  two  signals  is  in  the  ratio  of  the  S  to  initial  P  wave  amplitude 
levels.  The  difference  here  is  also  roughly  a  factor  of  three  with  relatively 
larger  S  to  P  ratios  for  STERLING  HE  than  for  STERLING.  Thus,  this 
observation  seems  to  be  consistent  between  station*  "OS  and  20S. 

In  an  attempt  to  develop  a  better  understanding  of  these  relative  P  and 
S  wave  differences  between  STERLING  and  STERLING  HE,  spectral 
analyses  were  performed  on  the  signals  recorded  at  stations  lOS  and  20S. 
Figure  3  shows  the  initial  P  wave  signal  and  pre-signal  noise  spectra 
determined  from  the  station  lOS  recordings  of  STERLING  (top)  and 
STERLING  HE  (bottom).  It  can  be  seen  from  these  figures  that  the  P  wave 
spectral  amplitude  levels  are  well  above  the  noise  levels  out  to  beyond  40  Hz 
in  both  cases.  Moreover,  the  signal  spectra  for  the  two  events  appear  to  be 
quite  comparable  with  both  showing  peaks  at  about  13  Hz,  above  which  they 
fall  off  slowly  toward  higher  frequencies.  The  corresponding  P  wave  signal 
and  noise  spectra  derived  from  the  station  20S  recordings  of  the  two  events 
are  ^own  in  Figure  4  where  it  can  be  seen  that,  while  the  signal-to-noise 
ratio  for  STERLING  is  again  good  out  to  about  40  Hz,  the  STERLING  HE 
signal  spectrum  is  above  noise  only  out  to  about  25  Hz.  In  particular,  a 
rather  prominent  noise  spike  is  apparent  in  the  STERLING  HE  spectra  at 
about  30  Hz  for  both  signal  and  noise.  Once  again,  the  P  wave  spectral 
shqres  for  the  two  events  appear  to  be  quite  similar  over  the  frequency 
range  in  which  they  can  be  validly  compared.  These  spectra  have  somewhat 
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Hgae  4.  Relative  P-wave  amplitude  spectra  at  station  20S  for  STERLING  (top)  and 
STERLING  HE  (bottom). 
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lower  comer  frequencies  than  the  corresponding  station  lOS  spectra  of 
Figure  3,  presumably  due  to  increased  effects  of  anelastic  attenuation  on  the 
propagating  P  wave  signals.  As  a  result,  the  P  wave  spectra  at  station  20S 
have  their  maximum  values  at  frequencies  near  4  to  5  Hz. 

The  corresponding  S  wave  spectra  for  STERLING  and  STERLING 
HE  are  plotted  in  Bgure  5  fm*  the  signals  at  station  lOS  (top)  and  at  station 
(bottom).  At  station  lOS  die  signal  spectrum  is  above  the  noise  over  the 
entire  frequency  band  out  to  50  Hz  for  both  STERLING  and  STERLING  HE 
although  the  STERLING  HE  signal  appears  to  be  close  to  die  noise  level  at 
frequencies  above  about  25  Hz.  At  station  20S  the  STERLING  S-wave 
signal  spectrum  is  above  noise  over  the  entire  frequency  range,  but  the 
STERLING  HE  signal  drops  to  the  noise  level  near  25  Hz.  For  the 
STERLING  HE  spectrum  the  noise  spike  at  30  Hz,  observed  in  the  P>wave 
results,  is  again  apparent  in  the  S*wave  spectrum.  At  lOS  the  S-wave  spectra 
for  STERLING  and  STERLING  HE  exhibit  a  broad  plateau  at  low 
frequencies  with  a  comer  near  14  Hz;  maximum  spectral  amplitudes  are  at  9 
Hz  and  4  Hz,  respectively.  Comparison  of  the  STERLING  and  STERLING 
HE  S-wave  spectra  at  station  lOS  reveals  that  the  STERLING  signals  are 
relatively  richer  in  high  frequencies  above  about  5  Hz.  At  20S  the  S-wave 
spectra  show  similar  behavior.  However,  at  20S  the  low-frequency  spectral 
plateau  seen  at  lOS  is  not  apparent;  this  can  again  probably  be  attributed  to 
de-emphasis  of  the  comer  by  higher  attenuation  at  the  more  distant  station. 
As  a  result,  the  S-wave  spectra  are  observed  to  decay  quite  rapidly  at 
frequencies  above  the  spectral  peak  which  is  near  5  Hz  for  STERLING  and 
near  4  Hz  for  STERLING  HE.  Again,  the  spectral  decay  is  less  rapid  for 
STERLING  so  that  the  S-wave  spectra  at  20S  are  enriched  at  high 
frequencies  relative  to  STERLING  HE.  Therefore,  at  both  stations  lOS  and 
20S  die  decoupled  STERLING  explosion  produced  S-waves  relatively  richer 
in  high  frequencies  than  those  from  the  tamped  chemical  explosion, 
STERLING  HE. 
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The  frequency  dependence  of  the  S/P  wave  ratios  for  the  two  source 
types  have  also  been  compared  through  narrowband  filter  analyses  of  the  P 
and  S  wave  signals  recorded  at  stations  lOS  and  20S.  In  this  approach,  the 
ratio  of  the  maximum  amplitudes  in  the  P  and  S  time  windows  for  each  filter 
output  are  computed  and  the  results  are  plotted  versus  the  center  frequency 
of  the  filter.  Such  ratios  are  plotted  in  Figure  6  for  station  lOS  (left)  and 
20S  (right),  where  the  results  for  STERLING  are  shown  as  solid  lines  and 
those  for  STERLING  HE  as  the  dashed  lines.  It  can  be  seen  that  in  all  cases 
the  S/P  ratios  decrease  with  increasing  frequency  from  maximum  values  of 
3-10  at  low  frequencies  to  values  on  foe  order  of  0.5-1.0  at  high  frequencies. 
These  results  corroborate  the  previously  noted  general  observations  made 
from  the  time  histories  and  individual  spectra  to  the  effect  that  the  S  wave 
signals  are  of  lower  dominant  frequency  than  foe  corresponding  P  wave 
signals  at  both  stations.  Of  greater  potential  significance  for  purposes  of 
event  identification,  however,  is  the  observation  that  the  S/P  ratios  are 
generally  about  a  factor  of  two  larger  for  STERLING  HE  than  for 
STERLING  over  rather  broad  frequency  bands  at  both  stations.  At  station 
lOS  this  difference  persists  over  the  entire  frequency  band  out  to  more  than 
20  Hz,  while  at  20S  the  two  ratios  are  different  out  to  about  12  Hz,  above 
which  they  appear  to  converge.  However,  this  apparent  high  frequency 
convergence  may  just  reflect  the  reduced  signal-to-noise  ratios  at  higher 
frequencies  at  this  more  distant  station.  In  any  case,  this  large,  frequency 
independent  difference  in  the  S/P  spectral  ratios  is  generally  consistent  at  the 
two  stations  and  its  persistence  over  such  a  broad  range  of  frequencies  is  a 
major  puzzle.  If  it  could  be  established  that  such  differences  are 
characteristic  of  CE  and  nuclear  cavity  decoupled  explosion  sources,  then  it 
would  be  possible  to  define  a  valuable  new  discriminant  for  use  in  seismic 
monitoring.  However,  before  proceeding  to  such  a  definition,  it  is  first 
necessary  to  establish  the  physical  plausibility  and  general  applicability  of 
any  such  discriminant.  Therefore,  it  is  appropriate  to  consider  the  source 
characteristics  and  near-source  environments  of  the  STERLING  and 
STERLING  HE  events  in  greater  detail. 
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As  was  noted  previously,  STERLING  and  STERLING  HE  were 
detonated  at  nearly  the  same  depth  in  the  Tatum  salt  dome  with  a  horizontal 
separation  of  about  350  m.  Their  relative  locations  with  respect  to  the 
boundary  of  tl^  salt  dome  at  source  depth  are  shown  in  Figure  7  where  it 
can  be  seen  that,  while  their  separation  is  not  negligible  with  respect  to  the 
geometry  of  the  boundary,  both  shot  points  are  located  more  than  400  m 
from  the  nearest  boundary  surface.  STERLING  HE  consisted  of  2.25  tons 
of  nitromethane  which,  according  to  Springer  et  al  (1968),  is  the  energy 
equivalent  of  about  2.7  tons  of  TNT.  The  explosive  was  distributed  in  a 
circular  cylinder  with  a  diameter  of  0.99  m  and  a  length  of  2.44  m. 
STERLING,  on  the  odier  hand,  was  a  380  ton  nuclear  explosion  in  a  roughly 
spherical  cavity  with  a  radius  of  about  17  m.  These  two  source  geometries 
are  shown  in  Figure  8  where  they  are  compared  with  that  of  the  recent 
"chemical  kiloton"  CE  test  conducted  at  NTS  in  September,  1993.  It  can  be 
seen  diat  none  of  these  initial  charge  geometries  are  spherically  symmetric 
and,  therefore,  the  question  arises  whether  such  variations  in  initial  source 
shape  could  be  responsible  for  the  observed  differences  in  S/P  ratio  between 
STERLING  and  STERLING  HE.  The  most  direct  approach  to  this  problem 
would  be  to  carry  out  a  nonlinear  finite  difference  simulation  of  the 
evolution  of  the  STERLING  HE  seismic  source  and  then  use  the  results  to 
quantitatively  evaluate  the  relative  efficiencies  of  S  and  P  wave  generation 
by  such  a  source.  However,  such  a  detailed  simulation  analysis  lies  outside 
the  scope  of  the  current  project.  Moreover,  there  is  much  past  experience 
which  suggests  that  such  a  simulation  would  reveal  only  very  minor 
departures  from  spherical  symmetry  in  the  resulting  seismic  source  function. 
That  is,  it  has  been  found  that  the  strong  nonlinear  response  of  the  medium 
surrounding  the  explosion  tends  to  compensate  for  any  asymmetries  in  the 
initial  charge  configuration,  producing  a  much  more  nearly  spherically 
symmetric  seismic  source  function.  Thus,  for  example.  Rimer  et  al.  (1993) 
in  their  nonlinear  simulation  of  the  chemical  kiloton  charge  geometry  shown 
in  Figure  8  found  that  the  simulated  final  cavity  volume,  cavity  pressure  and 
yield  radius  were  virtually  identical  to  the  corresponding  values  computed 
for  a  spherically  symmetric  charge  distribution  having  the  same  volume. 


14 


t 


K 


500  a 


Rguie7.  Source  locations  of  STERLING  and  STERLING  HE  with  respect  to  the 
boundary  of  the  Tatum  salt  <k)me  at  source  depth. 
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Similar  results  were  obtained  by  Stevens  et  al  (1993)  in  their  simulation 
analyses  of  the  seismic  source  functions  corresponding  to  some  strongly 
asymmetric  nuclear  charge  configurations  and  by  Stevens  et  al.  (1991)  in 
their  simulations  of  the  seismic  respond  of  ellipsoidal  cavities  having  much 
larger  aspect  ratios  than  STERLING  HE.  Therefore,  we  conclude  that  it  is 
highly  unlikely  that  the  initial  charge  asymmetry  of  STERLING  HE  is 
responsible  for  the  observed  differences  in  S/P  ratios  with  respect  to 
STERLING. 

An  alternate  possibility  is  that  differences  in  propagation  path 
characteristics  might  be  responsible  for  the  observed  variations  in  S/P  ratios. 
Thus,  although  STERLING  and  STERLING  HE  were  detonated  within  350 
m  of  one  another  in  the  same  salt  dome,  their  positions  relative  to  the  dome 
boundary  were  somewhat  different,  as  illustrated  in  Figure  7.  Thus,  it  is 
possible  that  the  interaction  of  the  outgoing  compressional  waves  with  the 
dome  boundaries  was  different  in  the  two  cases,  resulting  in  differences  in 
the  P  to  S  conversion  efficiency.  A  remaining  difficulty  with  this  hypothesis 
is  that  it  is  not  dear  why  any  such  boundary  effect  should  be  nearly  constant 
over  such  a  broad  frequency  band.  These  difficult  issues  will  ultimately 
have  to  be  addressed  using  three-dimensional  finite  difference  simulations  of 
the  two  explosions.  In  any  case,  it  is  probably  instructive  to  note  the  rather 
significant  seismic  differences  which  can  be  observed  at  a  common  station 
from  similar  sources  in  such  close  proximity  to  one  another. 

In  summary,  near-regional,  broadband  seismic  data  recorded  from  the 
STERLING  nuclear  cavity  decoupling  test  and  the  nearby  STERLING  HE 
tamped  CE  test  have  been  carefully  compared  and  evaluated  for  differences 
which  might  be  diagnostic  of  source  type.  It  has  been  demonstrated  that 
although  the  relative  spectral  compositions  of  the  P  and  S  waves  observed 
from  the  two  sources  appear  to  be  roughly  comparable,  the  S/P  spectral 
ratios  are  significantly  different  over  a  broad  frequency  band  extending 
from  at  least  1  to  10  Hz.  Evaluation  of  the  possible  effects  of  differences  in 
initial  charge  geometry  have  been  evaluated  in  a  preliminary  fashion  and  it 
has  been  concluded  that  it  is  unlikely  that  such  effects  can  explain  the 
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observed  differences  in  S/P  ratios  in  this  case.  An  examination  of  the 
potential  Influence  of  variations  in  the  locations  of  the  two  explosions  with 
req)ect  to  the  boundaries  of  the  Tatum  salt  dome  suggests  that  this  might  be 
a  contributing  factor,  although  it  is  not  clear  at  this  time  that  such  a 
mechanism  can  explain  the  frequency  independent  nature  of  the  observed 
seismic  differences.  In  conclusion,  additional  experimental  data  will  have  to 
be  evaluated  before  it  can  be  determined  whether  observed  seismic 
differences  of  the  kinds  noted  for  STERLING  and  STERLING  HE  can  be 
used  to  reliably  discriminate  between  tamped  CE  and  nuclear  cavity 
decoupled  explosions. 


3.  ComiMrisoii  of  P  Wave  Signals  Recorded  at  NORSAR  from 
CE  and  Nuclear  Cavity  Explosions  in  the  North  Caspian 
Basin 


The  area  surrounding  the  Caspian  Sea  provides  a  good  test  bed  for 
evaluating  procedures  for  identif^ng  cavity  decoupled  nuclear  explosions  in 
that  it  contains  the  Azgir  test  site  at  which  the  Russians  have  acknowledged 
conducting  a  number  of  nuclear  explosions  in  salt  cavities.  Moreover, 
although  the  natural  seismicity  of  the  area  surrounding  the  Azgir  site 
iqf^ars  to  be  quite  low,  there  appear  to  be  significant  numbers  of  CE  events 
occurring  in  the  salt  dome  regions  of  the  Nordi  Caspian  Basin  which  can 
serve  as  a  source  of  seismic  data  to  be  compared  with  corresponding  data 
from  nuclear  cavity  explosions  fra:  evaluation  of  potential  discriminants.  In 
the  following  discussion,  we  will  first  briefly  review  the  seismicity  of  die 
Caspian  region  and  then  proceed  to  a  comparison  of  NORSAR  P  wave  data 
recorded  from  selected  CE  and  nuclear  cavity  explosions  in  the  North 
Caspian  Basin. 

Strong  differences  in  tectonics  and  natural  seismicity  exist  in  the 
region  surrounding  the  Caspian  Sea.  The  southern  region,  including  the 
Caucasus  mountains  west  of  the  Caspian  and  die  Kopet  Dagh  area  east  of  the 
Caspian,  is  tectonically  active  with  a  high  level  of  natural  earthquake 
activity.  In  contrast  the  northern  Caspian  basin  is  a  stable  platform  region 
with  little,  if  any,  natural  seismicity.  This  dichotomy  is  illustrated  in  Figure 
9  which  depicts  epicenters  of  earthquakes  with  energy  class  K  >  11  (mb  >  4) 
reported  in  d^  annual  Russian  earthquake  catalogs  for  the  period  1962  to 
1979.  The  earthquake  activity  tends  to  be  clustered,  but  the  map  shows  a 
clear  diminution  in  the  level  of  seismicity  north  of  about  44°  latitude  on  the 
west  side  of  the  Caspian  and  a  similar  drop  off  in  activity  north  of  40° 
latitude  on  the  east  side  of  the  Caspian.  Although  this  northward  drop  in 
seismicity  may  to  some  extent  reflect  Russian  reporting  practice  and  bias 
associated  with  the  locations  of  Soviet  seismograph  stations,  the  stable 
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geologic  conditions  and  total  absence  of  earthquakes  in  the  northern  Caspian 
basin  as  reported  in  international  catalogs  appears  to  confirm  that  the  area  is 
seismically  quiet.  This  is  potentially  significant  in  the  current  context,  since 
the  principal  Soviet  nuclear  testing  area  in  salt  was  located  in  the  northern 
Caspian  basin  at  Azgir  (cf.  Figure  9).  Thus,  this  region  provides  a 
prototype  for  the  environment  in  which  seismic  monitoring  for  cavity 
decoupled  explosions  might  take  place. 

At  first  glance,  the  results  summarized  in  Figure  9  might  suggest  that 
detection  of  a  seismic  event  in  the  vicinity  of  Azgir  is  equivalent  to  its 
identification  as  a  nuclear  explosion.  However,  even  though  this  region 
seems  to  be  largely  devoid  of  natural  earthquake  activity,  a  number  of  low 
level  seismic  events  are  in  fact  detected  in  the  area  every  year.  For  example, 
Ringdal  (personal  communication,  1984)  has  identified  more  than  180  events 
in  the  NORSAR  catalog  which  were  located  in  the  aseismic  area  bounded  by 
4S°N  and  52°N  north  latitude  and  by  42°E  and  55°E  longitude  during  the 
time  period  October,  1971  to  October,  1981.  The  NORSAR  locations  for 
these  events  are  shown  in  Figure  10  where  it  can  be  seen  that  they 
completely  surround  the  Azgir  nuclear  testing  area.  Ringdal  has  concluded 
that  most  of  these  events  are  probably  chemical  explosions  on  the  basis  of 
their  strong  clustering  during  the  local  daylight  hours.  This  fact  is 
illustrated  in  Figure  11  which  shows  a  comparison  of  the  distribution  of 
times  of  occurrence  of  the  seismic  events  from  Ringdal's  north  Caspian  list 
with  that  for  natural  earthquakes  from  the  Caucasus  area  to  the  south.  It  can 
be  seen  that  the  north  Caspian  events  from  the  NORSAR  list  are  strongly 
clustered  in  the  time  period  from  about  6  to  18  hours  GMT,  in  contrast  to 
the  Caucasus  events  which  seem  to  be  evenly  spread  throughout  the  24  hour 
period.  This  strong  peaking  of  the  events  from  the  NORSAR  list  suggests 
that  the  majority  of  them  are  not  natural.  Given  this  observation  and  the 
presence  of  oil-shale  development  and  construction  activities  in  the  area,  it 
seems  reasonable  to  conclude  that  these  events  are  mostly  chemical 
explosions  similar  to  those  which  might  be  found  in  other  seismically  quiet 
areas  of  the  world. 
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Hgnie  1 1.  Cmnpariscm  of  the  distributitn  of  times  of  occurrence  of  the  seismic  events  from 
Ringdal's  nnth  Caqjian  list  (top)  with  that  for  natural  earthquakes  from  the  Caucasus  area 
to  the  south  (bottom). 
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As  was  iodicated  above,  a  number  of  North  Caspian  nuclear  and 
pfesumed  CE  events  have  been  detected  by  the  NORSAR  array  and  listed  in 
the  NORSAR  event  bulletin.  Table  1  lists  the  sample  of  12  such  events 
which  have  been  selected  for  comparative  analysis.  It  can  be  seen  that  the 
first  eight  of  these  have  been  designated  as  presumed  CE  events,  while  the 
last  four  have  been  identified  by  the  Russians  as  having  been  nuclear 
explosions  in  salt  cavities  at  Azgir  (Sultanov,  personal  communiciuion, 
1993).  It  can  be  seen  that  the  NORSAR  mb  estimates  listed  in  this  table 
indicate  that  the  samples  for  the  two  source  types  encompass  about  the  same 
range  in  magnitude,  thereby  permitting  direct  comparisons  of  the  recorded 
seismic  data.  The  NORSAR  locations  of  these  selected  events  are  plotted  on 
the  map  of  the  North  Caspian  region  in  Figure  12  where  it  can  be  seen  that 
the  estimated  locations  of  the  presumed  CE  events  surround  the  locations  of 
the  nuclear  explosions.  Note  from  Table  1  and  Figure  12  that  Events  #1  • 
#3  appear  to  be  co-located  and  were  detonated  on  the  same  day  at  intervals 
of  about  14  seconds,  which  would  seem  to  confirm  that  at  least  these  diree 
are  CE  events.  It  is  difficult  to  assess  the  accuracy  of  the  NORSAR  locations 
given  the  available  information,  but,  given  the  size  of  the  events  and  their 
distance  from  NORSAR  (A  -  25^),  it  would  not  be  surprising  if  the 
epicentral  uncertainties  are  on  the  order  of  100  km  or  more.  This  rough 
estimate  agrees  with  the  observed  maximum  mislocation  of  the  Azgir 
nuclear  explosions,  which  is  on  the  order  of  125  km  for  Event  #12. 

With  regard  to  the  nuclear  explosions,  the  Russians  have  confirmed 
that  explosions  #9  >  #12  in  Table  1  were  detonated  in  two  different  cavities 
in  salt  at  the  Azgir  test  site.  More  specifically,  Adushkin  et  al  (1992)  have 
indicated  that  Event  #9  was  an  8  kt  decoupled  nuclear  test  conducted  in  the 
38  m  radius,  air-filled  cavity  produced  by  the  64  kt  tamped  nuclear 
explosion  conducted  at  Azgir  on  12/22/71.  Similarly,  Sultanov  et  al 
(personal  communication,  1993)  report  that  Events  #10  -  #12  were  sub- 
kiloton  nuclear  tests  conducted  in  the  32  m  radius,  water-filled  cavity 
produced  by  the  25  kt  tamped  nuclear  explosion  of  7/1/68.  Thus,  only  Event 
#9  was  a  true  decoupling  test  and  even  diat  was  only  partially  decoupled  as  a 
result  of  the  fact  that  its  yield  was  a  factor  of  three  or  more  larger  than  that 
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Table  1. 

NORSAR  Source  Parameters  of  Selected 
North  Caspian  Events 


Event#* 

Date 

Origin  TimedJD 

NUt 

E  Long. 

mb 

1(CE) 

09/14/78 

04:00:02 

51.2 

53.2 

3.6 

2(CE) 

09/14/78 

04:00:16 

51.2 

53.2 

3.6 

3(CE) 

09/14/78 

04:00:29 

51.2 

53.2 

3.6 

4(CE) 

04/28/84 

14:53:48 

46.2 

45.8 

4.5 

5(CE) 

06/27/86 

07:25:44 

48.1 

48.7 

4.1 

6(CE) 

07/25/86 

07:22:35 

47.4 

49.9 

3.9 

7(CE) 

08^/86 

06:31:37 

48.3 

49.4 

4.0 

8(CE) 

09/28/86 

14:08:16 

45.6 

47.5 

4.4 

9(NO 

03/29/76 

07:00:03 

48.0 

48.0 

3.9 

10  (NO 

10/14/77 

07K)0:02 

48.0 

48.0 

3.9 

11  (NO 

09/12/78 

05:00K)1 

48.0 

48.0 

3.5 

12  (NO 

01/10/79 

07:59:53 

46.9 

48.3 

4.1 

*  CE  deooies  presumed  chemical  explosion  and  NC  denotes  nuclear  cavity 
test 
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Hgm  12.  NORSAR  locations  of  the  selected  Caspian  events.  Open  circles  denote  nuclear 
caviQf  tests  ndttle  astetidcs  denote  (uesumed  chemical  explosions. 
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specified  by  the  Latter  criterion  for  full  decoupling  in  a  cavity  of  that  size. 
Mtxeovtt,  the  water-filled  cavity  tests  were  observed  to  have  coupled  better 
than  tamped  explosions  of  the  same  yield  in  salt,  consistent  with  previous 
experience  of  the  relative  seismic  coupling  efficiency  of  nuclear  explosions 
in  water.  However,  all  four  of  these  nuclear  tests  were  conducted  in  nearly 
spherical  cavities  in  salt  and,  therefore,  except  for  absolute  amplitude  level, 
they  might  be  expected  to  produce  low  frequency  seismic  signatures 
comparable  to  those  expected  from  fully  decoupled  nuclear  explosions. 

Full  NORSAR  array  data  have  been  retrieved  for  each  of  the  12 
events  of  Table  1  and  previewed  to  assess  relative  data  quality  across  the 
array.  In  agreement  with  the  results  of  previous  studies  by  Ringdal  et  al 
(1983),  it  was  found  that  subairay  02B  consistently  provides  the  highest 
quality  signals  frxr  diese  north  Caspian  events  and,  consequently,  our  analyses 
have  focused  on  data  recorded  at  the  six  elements  of  this  subarray.  The 
recordings  at  the  center  element  (i.e.,  O2B0)  for  the  12  selected  events  are 
compared  in  Figure  13.  It  can  be  seen  that  the  signal-to-noise  ratios  at  this 
station  are  generally  quite  good  for  these  events,  although  there  is  some 
evident  interference  from  longer  period  microseismic  noise.  Due  to  the 
differences  in  relative  frequency  content,  this  noise  can  be  effectively 
suppressed  by  filtering  and  Hgure  14  shows  the  results  of  filtering  the 
seismograms  of  Figure  13  through  a  bandpass  filter  with  comer  frequencies 
at  0.8  and  6.0  Hz.  It  is  evident  from  these  figures  that  there  are  some 
significant  differences  in  the  short-period  P  waves  produced  by  these  events. 
Fot  exaiiq>le,  the  recordings  from  the  HE  Events  #4,  #5  and  #6  are  more 
emer^nt  and  have  higher  coda  levels  than  those  of  the  nuclear  cavity  Events 
#9  -  #12,  which  are  very  similar  to  one  another.  On  the  other  hand,  the 
rec(H*dings  from  some  of  the  HE  Events  (e.g.,  #1  -  #3)  appear  to  be  quite 
similar  in  character  to  the  nuclear  recordings.  Thus,  in  terms  of  general 
time  domain  characteristics,  the  recordings  from  the  presumed  chemical 
explosions  show  a  degree  of  variability  which  appears  to  encompass  the 
characteristics  of  the  corresponding  decoupled  seismograms. 
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#7 

»8 

«9 


«10 

*11 

*12 


5  seconds 
I - 1 


FifvelS.  NORSARsttbairaycleincntOZBO  recordings  of  selected  Caspian  seisi^ 

eveitfs. 
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5  seconds 


Rgoie  14.  Hllered  (0.8-6.0  Hz)  NORSAR  subanay  element  O2B0  recordings  of  selected 
Q^nan  events. 
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The  next  step  in  the  analysis  proceduie  was  to  estimate  P  wave  signal 
and  noise  spectra  for  the  12  events  at  each  of  the  six  elements  of  subarray 
02B.  The  data  were  recorded  at  a  sampling  rate  of  20  samples/second  and 
the  spectra  were  computed  using  the  first  4.5  seconds  of  the  P  wave  signal 
and  a  corresponding  window  of  pre-signal  noise.  The  resulting  spectra  were 
smoothed  and  cmrected  for  instrument  response  and  representative  spectra 
for  Events  #1,  #4  and  #10  are  displayed  in  Figures  15,  16,  and  17, 
respectively.  With  reference  to  Figures  15-17,  it  can  be  seen  that  the  usable 
passband  generally  extends  from  about  1-2  Hz  to  about  7-9  Hz.  Although 
there  is  some  evidence  of  site  effects  on  the  spectra  observed  at  the 
individual  elements,  the  P  wave  spectral  composition  has  been  found  to  be 
quite  stable  across  the  02B  subarray.  This  is  illustrated  in  Figures  18-20 
which  show  the  individual  element  and  subarray-averaged  P  wave  spectra 
corresponding  to  the  events  of  Figures  15-17,  plotted  over  the  passbands 
containing  usable  signal-to-noise  ratios.  The  dashed  lines  enveloping  the 
subarray  averages  in  these  figures  denote  the  one  standard  error  of  estimate 
levels  which,  for  the  six  element  array,  correspond  roughly  to  95  percent 
confidence  intervals  about  the  average  values.  For  the  HE  events 
exenqilified  by  Figures  18  and  19,  the  standard  error  of  estimate  over  this 
frequency  band  averages  to  less  than  a  factor  of  1.5.  Somewhat 
surprisingly,  the  variation  in  P  wave  spectra  across  the  subarray  appears  to 
be  noticeably  larger  for  the  nuclear  cavity  explosions,  averaging  to  nearly  a 
factor  of  2.0,  as  illustrated  in  Figure  20.  There  is  no  obvious  reason  for  this 
difference  and  it  may  just  reflect  the  fact  that  the  decoupled  events  are 
generally  older  and  may  have  been  recorded  during  a  period  in  which  the 
instrument  calibration  factors  were  less  well  determined  than  they  have  been 
in  more  recent  times. 

Examination  of  these  subarray-averaged  P  wave  spectra  indicates  that 
the  events  can  be  roughly  divided  into  three  groups,  consisting  of  Events  #1 
-  #3,  Events  #4  -  #8  and  the  nuclear  cavity  Events  #9  -  #12.  The  subarray- 
averaged  spectra  within  each  of  these  three  groups  are  conq)ared  in  Figures 
21-23,  respectively,  where  the  individual  spectra  have  been  normalized  to 
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Ralativ*  Spectral  Aaplituda 


Selativ*  Spectral  Aaplitiide 


ia3<»67i»10  13I4S6789  10 

Ht  t.  Hr 


Hgae  16.  NCHISAR  subainy  02B  P  wave  signal  (solid  line)  and  noise  (dashed  line) 
q)ectrif(tf  Event#4. 
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telative  Spectral  Aaplitude 


Rgnie  17.  NORSAR  subairay  02B  P  wave  signal  (solid  line)  and  noise  (dashed  line) 
qiectia  for  Event  #10. 
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Rela 


iteire  22.  OM^MBtisoa  of  NCH(SAR  02B  sobamy-avenged  P  wave  spectia  for  Events 
#4-#8. 
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Spectral  Anplitude 


Rgiiie23.  Compiriscm  of  NCXlSAR02Bsiibiiny-avenigedP  wave  spectn  for  Events 


IIk  tame  average  amplitude  level  over  the  available  passband.  It  can  be  seen 
that  the  spectral  shapes  within  each  group  are  quite  consistent  over  the 
frequency  band  available  for  analysis.  Thus,  for  example,  the  spectra  for 
nuctear  Events  #10  and  #1 1  in  Figure  23  are  essentially  identical,  as  might 
be  expected  since  they  are  of  similar  size  and  were  detonated  in  die  same 
cavity.  More  surprising  is  the  close  similarity  of  the  P  wave  spectral  shi^s 
for  die  HE  Events  #4  •  #8  in  Figure  22,  particularly  considering  the  fact  that 
diey  leptesent  a  range  of  source  size  and  testing  location.  Hgure  24  shows  a 
comparison  of  the  average  spectra  for  each  of  the  three  groups  of  Figures 
21-23,  approximately  normalized  to  the  same  low  frequency  level.  It  can  be 
seen  diat  die  group  consisting  of  the  [Resumed  chemical  explosions  #4  -  #8 
rqaiesents  the  lowest  frequency  source,  while  the  group  consisting  of  the 
presumed  chemical  explosions  #1  -  #3  iqipears  to  represent  the  broadest  band 
source.  Once  again,  as  with  the  time  domain  characteristics,  die  average  P 
wave  spectra  of  the  presumed  chemical  explosions  in  die  two  groups  differ 
mme  from  <me  another  than  eitho*  does  from  the  average  nuclear  cavity  P 
wave  spectrum.  This  observation  is  illustrated  more  specifically  in  Figure 
25  where  the  ratios  of  the  average  spectra  for  the  two  HE  groups  of  events 
computed  with  respect  to  spectrum  for  cavity  decoupled  Event  #9  are 
compared.  This  figure  concisely  summarizes  the  principal  conclusion  of  the 
NORSAR  analysis,  which  is  that  the  variation  of  the  characteristics  of  the 
recorded  P  waves  from  presumed  chemical  explosions  envelopes  the 
characteristics  of  the  recorded  P  waves  from  nearby  nuclear  cavity 
explosions  of  conqiarable  magnitude.  Thus,  this  evidence  suggests  that 
teleseismic  P  wave  data  in  the  band  frtim  about  1  to  10  Hz  is  unlikely  to 
provide  a  reliable  basis  for  discriminating  between  these  two  source  types. 
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Rela 


£,  Hz 


Rgue  24.  Comparison  of  NORSAR  02B  P  wave  spectral  average  for  Events  #l-#3  with 
corrmponding  averages  fOT  Events  #4-#8  and  Events  t^l2. 
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4.  Comparison  of  Broadband  Near-Regional  Seismic  Record 
ii^  from  Soviet  CE  and  Cavity  Decoupled  Explosions  in 
the  Caspian  Sea  Region 


As  was  noted  in  Section  3,  we  currently  have  knowledge  of  only  two 
full  scale  nuclear  decoupling  experiments;  namely,  the  U.S.  test  STERLING 
and  the  Soviet  air-filled  cavity  test  at  Azgir.  The  Soviet  cavity  decoupling 
experiment  was  carried  out  on  3/29/76  when  an  8  kt  nuclear  device  was 
detonated  at  a  depth  of  987  m  in  a  38  m  radius  cavity  in  salt  produced  by  a 
previous  64  kt  tan^wd  nuclear  explosion  conducted  at  that  site  on  12/22/71 
(Adushkin  et  al.^  1992).  Broadbami  seismic  data  were  recorded  from  this 
decoupled  test  at  near-regional  distances  ranging  from  about  1  to  113  km 
and  these  analog  data  were  subsequently  digitized  at  a  rate  of  200 
sanq>les/second,  converted  to  ground  motion  amplitudes  and  installed  as  a 
database  at  the  ARPA  CSS  (personal  communication,  I.  Kitov,  1992).  These 
data  provide  an  important  reference  base  for  conq>arison  with  seismic  data 
recorded  in  the  same  distance  range  from  CE  events  in  the  Caspian  Sea 
region.  However,  before  conq)arisons  can  be  made  for  the  nominal,  fiilly 
decoupled  1  kt  evasion  scenario,  it  is  necessary  to  frrst  normalize  the 
obs^ed  data  to  account  for  the  fact  that  the  Azgir  decoupling  experiment 
was  at  a  yield  of  8  kt  and  was  only  partially  decoupled.  This  has  been 
accoiiq>lished  using  the  approximate,  theoretical  scaling  model  described 
below. 

Let  Zi(t),  Z2(t)  denote  the  seismograms  for  a  particular  component  of 
motion  recorded  at  a  fixed  station  from  explosions  at  the  same  source 
location,  with  corresponding  Fourier  transforms  Zi(®),  Z2(®).  Then,  in 
operational  form,  the  ground  motion  spectra  can  be  written  as 
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Zl(OJ)  a  Sl(<0)  T((D) 


(1) 


Z2«D)  =  S2((0)  T((D) 


where  Si(co),  S2(<ii>)  denote  the  seismic  source  functions  corresponding  to  the 
two  «q>losions  and  T(<o)  is  the  propagation  path  transfer  function  which  is 
common  to  both.  It  dien  follows  diat,  to  this  degree  of  ^proximation 


Z2((D)  = 


S2(co) 

SiCca) 


Zl((D) 


(2) 


oo 

Z2(t)  =  J  Zl(x)  S(t  - 1)  dt  (3) 

•«o 


where  S(t)  is  the  explosion  source  scaling  operator.  Thus,  the  expected 
ground  motion  time  history  for  the  second  explosion  can  be  determined 
from  the  observed  ground  motion  from  the  first  explosion  if  the  spectral 
noio  of  the  two  seismic  source  functions  can  be  estimated.  For  the  simplest, 
si^mically  synunetric  model  of  the  explosion  source,  the  required  source 
spectral  ratio  has  the  form  (Mueller  and  Murphy,  1971) 

^  _  P2(a.)reb  «»°.  ^  ^ 

Si(©)  Pl(<h)  reli  10^2  +  i  ©o2®  - 


where  reli,  rel2  are  the  elastic  radii  of  the  two  sources  at  which  spherically 
symmetric  pressures  pi(o>),  p2(<o)  act  and 


a 


(5) 


X+2\l 

4^l 


with  a  die  P  wave  velocity  and  X,  ^  the  Lame  constants  characteristic  of  the 
common  source  medium.  Now  for  explosions  at  the  same  depth  in  a  fixed 
medium,  the  elastic  transition  pressure  should  be  constant  then,  assuming  a 
step  function  i^proximation  for  the  pressure  profiles  acting  at  the  elastic 
radii,  it  follows  that  p2((o)  =  pl((*>)  in  (4)  and  the  modulus  of  the  source 
spectral  ratio  can  be  written  simply  as 


|S«D)|  = 


82(0)  rel2 
Si((0)  *  reli 


. . 

(CObl  -  pto2)2  +  (0oj(02 
(£062  • 


(6) 


Thus,  the  !sw  and  high  frequency  asymptotic  values  of  the  source  spectral 
ratio  are  given  by 


lim 


|S(©)|  = 


,  2 
rcl2  ©ol 

r 

reli  ©02 


lim 

09  oo 


|S(©)|  = 


rel2 

reli 


(7) 


OT,  since  rel  is  jvoportional  to  the  cube  root  of  the  yield  W  for  explosions  at 
a  fixed  depth  in  a  given  medium. 
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(8) 


This  is  the  familiar  cube  root  scaling  model  for  the  explosion  seismic  source 
function. 


The  time  domain  operator  corresponding  to  Equation  (6)  for  the  case 
X  =  p  has  the  form 

S(t)  =  -^e  3  ^Cicos^^t0o2t  +  C2sin  ^^tDo2tj  +  ^8(t) 

(9) 

where 


Then  ^fining  S(t)  =  Si(t)  +  (rel2/reli)  5(t)  in  (9),  we  can  write 


OO 

Z2(0  =  j  Zi(T)  S|(t-t)  dT  +  ^  Zi(t)  (11) 

-OO 

It  follows  that,  given  an  observed  seismogram  from  an  explosion  with 
known  source  characteristics,  the  seismogram  to  be  expected  at  that  same 
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station  from  an  explosion  with  different  source  characteristics  in  tl^  same 
region  can  be  approximated  using  scaling  relations  of  the  form  exemplified 
by  Equations  (9)  -  (11). 

The  sinq)le  scaling  procedures  described  above  have  proven  to  work 
quite  well  in  a  number  of  different  applications.  For  example.  Figure  26 
shows  the  results  of  scaling  observed  STERLING  data  from  near-regional 
station  lOS  at  a  range  of  16  km  to  the  source  conditions  of  SALMON  using 
the  source  parameters  inferred  by  Murphy  (1969)  for  these  two  explosions. 
It  can  be  seen  that,  despite  the  sinq>licity  of  the  source  approximations,  the 
agreement  between  the  scaled  STERLING  and  observed  SALMON  seismic 
data  at  this  station  is  quite  remarkable.  This  would  seem  to  confirm  the 
validity  of  die  procedure  as  applied  to  co-located  events  in  the  same  medium. 
A  remaining  question  concerns  the  applicability  of  SALMON/STERLING 
seismic  source  scaling  and  decoupling  relations  to  explosions  in  the  Azgir 
salt  domes.  Although  additional  work  is  needed  to  fully  assess  the 
consistency  of  the  seismic  source  coupling  in  these  different  salt 
environments,  some  recently  acquired  Soviet  free-field  data  from  a  tamped 
explosion  at  Azgir  suggests  that  they  might  be  quite  similar.  These  new 
free-field  data  were  recorded  from  the  1.1  kt  tamped  nuclear  explosion 
which  was  conducted  at  a  depth  of  161  m  at  Azgir  on  4/22/66.  For  this  test, 
vertical  component  particle  velocity  waveforms  were  measured  above  the 
detonation  point  at  ranges  extending  from  63  to  161  m  in  a  satellite  hole 
offset  from  the  emplacement  hole  by  about  10  m  (personal  communication, 
I.  Kitov,  1993).  Using  simple  cube-root  scaling  as  a  first  approximation,  the 
receding  made  at  a  range  of  84  m  above  shot  depth  corresponds  to  a  range 
of  142  m  at  the  5.3  kt  yield  of  SALMON.  The  cube-root  scaled  version  of 
this  recording  is  shown  in  Figure  27  where  it  is  compared  with  the 
SALMON  shot  depth  recording  at  a  range  of  166  m.  It  can  be  seen  that  the 
arrival  times,  waveshapes  and  amplitude  levels  (preliminary  for  Azgir)  are 
remarkably  similar  for  these  two  recordings,  which  suggests  that  the  seismic 
source  coupling  characteristics  of  Azgir  salt  may  not  be  very  different  from 
those  of  SALMON  salt.  Additional  data  will  have  to  be  analyzed  in  order  to 
definitively  test  this  hypothesis.  However,  these  preliminary  data  suggest 
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26.  Comparison  of  observed  SALMON  and  scaled  STERLING  (SALMCH4  synthetic)  radial  compmieiit 
don  lOS. 
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Figufe  27.  Cbnqnrison  of  observed  SALMON  free-field  particle  velocity  recording  at  a 
range  of  166  m  widi  the  cube-root  scaled  Azgir  4/22/66  ftee-field  particle  velocity  recording 
at  a  scaled  range  ( at  5.3  kt )  of  142  m. 
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that  the  SALMON/STERLING  source  scaling  procedures,  illustrated  by  the 
results  shown  in  Figure  26,  should  provide  a  reasonable  first  {^{n-oximation 
f(v  ^>plication  to  Azgir  explosion  data. 

The  vertical  consonant  displ^ment  waveforms  recorded  from  the 
Azgir  cavity  decoupling  test  of  3/29/76  in  the  near-regional  distance  range 
extending  from  about  18  to  113  km  are  shown  in  Figure  28.  It  is  evident 
that  these  data  are  quite  complex  and  variable  as  a  function  of  distance. 
Thus,  for  example,  it  can  be  seen  that  while  the  first  arrival  at  2S.S  km  is 
inq>ulsive  and  simple,  that  at  the  nearby  23.0  km  station  is  much  more 
emergent  and  complex.  Similar  comments  apply  to  the  observations  at 
distances  greater  than  40  km,  which  show  no  obvious  explosion-like 
characteristics  of  the  kind  which  frequently  distinguish  teleseismic  P  waves 
recorded  from  underground  explosion  sources.  This  observed  high  degree 
of  variability  is  typical  of  broadband,  near-regional  recordings  from 
earthquakes  and  explosions  and  contributes  to  the  difficulty  of  defining 
reliable  discriminants  for  application  in  this  distance  range. 

Now,  before  the  data  of  Figure  28  can  be  scaled  to  that  expected  for 
the  1  kt  fully  decoupled  scenario,  it  is  necessary  to  first  define  an 
approximate  seismic  source  model  for  this  8  kt  partially  decoupled 
explosion.  Detailed,  nonlinear  finite  difference  simulations  of  this  source 
are  currently  underway  but,  for  the  present  application,  we  will  employ  the 
simplest  analytic  approximation  corresponding  to  a  step  in  pressure  acting  at 
the  elastic  radius.  Preliminary  modeling  of  the  near-field  P  wave  spectra 
for  this  test  suggests  an  equivalent  elastic  radius  of  about  50  m,  a  value 
which  is  significantly  larger  than  the  initial  cavity  radius  of  38  m.  This  is 
not  surprising,  since  a  cavity  with  a  radius  of  more  than  45  m  would  be 
required  to  match  STERLING  conditions  for  8  kt  at  this  depth.  It  follows 
that  cube-root  scaling  to  the  comparable  conditions  of  1  kt  in  a  19  m  radius 
cavity  at  die  depth  of  the  Azgir  cavity  should  give  an  elastic  radius  of  about 
25  m.  The  source  spectral  ratio  corresponding  to  Equation  (6)  for  this  8  kt 
to  1  kt  scaling  operator  is  shown  as  a  dashed  line  in  Figure  29,  where  it  can 
be  seen  diat  the  low  frequency  level  is  equal  to  0.125  (i.e.,  (25/50)^)  and  the 
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17.8  km  \  r\/vlA/VlAA/y/'nA/*v\^^^  6.28x10'^ cm 


Figure  28.  Wrdcal  component  ground  motionii  recorded  from  the  8  kt  Azgir  cavity  decoupling  teat  of  3/29f76. 


high  fluency  level  is  equal  to  0.5  (i.e.,  25/50)  in  accord  with  the 
asyiiq>totic  relations  given  by  Equations  (8).  Now  the  observed  (^coupling 
factor  for  the  3/29/76  Azgir  test  was  about  a  factor  of  15  (Adushkin  et  al., 
1992),  as  ofqiosed  to  the  observed  STERLING  nominal  "full"  decoupling 
factor  of  70.  Therefore,  in  order  to  approximately  scale  to  STERLING 
decoupling  conditions  for  1  kt  at  Azgir,  the  source  scaling  operator 
rqnesented  by  die  dashed  line  in  Figure  29  has  been  multiplied  by  the  ratio 
15/70  to  obtain  the  final  source  scaling  operator  shown  as  die  solid  line  in 
Figure  29.  This  is  the  operator  which  has  been  used  to  scale  the  near- 
regional  data  observed  from  the  8  kt  partially  decoupled  Azgir  test  to 
approximate  the  ground  motions  which  would  be  expected  from  a  1  kt  fully 
decoupled  explosion  at  that  same  source  location.  Note  that  for  frequencies 
less  than  about  15  Hz,  this  operator  in  nearly  constant,  independent  of 
frequency,  at  a  level  given  by  the  low  frequency  asymptotic  value  of 
Equaticm  (8).  It  follows  that  since  the  observed  8  kt  decoupled  Azgir  data  of 
Figure  28  are  dominated  by  components  having  frequencies  less  than  about 
10  Hz,  the  expected  effect  of  scaling  to  the  1  kt  fully  decoupled  scenario  is  to 
reduce  the  anqilitude  level  by  about  a  factor  of  40  without  significantly 
modifying  the  time  dependent  character  of  the  signals.  This  expectation  is 
confirmed  in  Figure  30  which  shows  the  approximations  to  the  1  kt  fully 
decoupled  ground  mcMions  obtained  by  ^plying  die  solid  line  source  scaling 
operator  of  Figure  29  to  the  observed  Azgir  waveforms  of  Figure  28. 
Conqiaiing  Figures  30  and  28,  it  can  be  seen  that  the  relative  spectral 
conqiosition  of  the  signals  is  unchanged  by  scaling  in  this  case,  in  contrast  to 
the  pronounced  change  for  the  SALMON/STERLING  scaling  results  shown 
previously  in  Figure  26.  The  inferred  peak  vertical  component  displacement 
amplitudes  from  Figure  30  are  plotted  versus  range  in  Figure  31  where  it 
can  be  seen  that  they  decrease  with  distance  at  an  average  rate  about  equal  to 
that  reported  by  Adushkin  et  al.  (1992)  for  tamped  explosions  at  Azgir. 

Unfortunately,  no  data  from  CE  events  have  been  recorded  at  the 
stations  used  to  record  the  Azgir  nuclear  decoupling  test.  However, 
numerous  mine  blasts  have  been  recorded  in  the  same  near-regional  distance 
range  at  IRIS  station  KIV  which  is  located  west  of  the  Caspian  Sea  at 
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Figure  30.  TheoietkaUy  scaled  veitkMl  componern  grouiKl  motkMis  conespooding  to  a  1  kt  fully 

decou|ded  exploskm  at  the  source  location  of  the  8  kt  Azgir  decoupling  lest  of  3/29/76. 


Kislovodsk,  approximately  600  km  southeast  of  Azgir.  The  location  of 
station  KIV  with  respect  to  Azgir  is  shown  in  Figure  32,  which  also  shows  in 
exfumded  scale  the  locations  of  two  mines  near  KIV,  as  well  as  the  locations 
of  five  satellite  stations  surrounding  KIV  which  have  recorded  seismic  data 
from  blasts  at  these  two  mines  (Rivi&re-Barbier,  1993).  These  stations 
recorded  data  at  an  extended  sampling  nue  of  125  samples/second  from  two 
well-documented  explosions  at  the  Tymyauz  mine  and  one  at  the  Zhako- 
Krasnogorskaya  mine  which  provide  a  broadband,  near-regional  seismic 
data  base  comparable  to  the  Azgir  nuclear  decoupling  data  recorded  in  the 
same  distance  range.  According  to  Rivi^re-Baibier  (1993),  the  two  well- 
documented  Tymyauz  mine  blasts  were  conducted  about  one  minute  i^art  on 
4/12/92.  Both  were  tipple-fired  explosions  detonated  in  multiple  boreholes 
extending  to  a  depth  of  15.6  m  below  the  surface  of  the  open  pit  mine.  The 
first,  denoted  4/12/92(1),  a  total  yield  of  27  tons  of  CE,  while  the  second, 
denoted  4/12/92(2),  had  a  total  yield  of  70  tons  CE.  The  documented  blast  at 
the  Zhako-Krasnogorskaya  mine  occurred  on  5/8/92  and  was  also  a  ripple- 
fired,  multiple  borehole  explosion  at  the  surface  of  this  gypsum  pit  mine, 
having  a  total  yield  of  1 1.9  tons  CE. 

Before  proceeding  to  waveform  comparisons,  it  is  appropriate  to 
examine  the  amplitude  levels  of  the  ground  motions  produced  by  these  mine 
blasts  as  they  relate  to  those  expected  in  the  same  distance  range  from  the 
nominal  1  kt  fully  decoupled  nuclear  explosion  at  Azgir.  The  peak  vertical 
displacements  observed  at  the  KIV  station  array  from  the  three  well- 
documented  mine  blasts  are  displayed  as  a  function  of  range  in  Figure  33. 
On  these  figures,  the  dashed  lines  represent  the  average  1  kt  fully  decoupled 
an^litude-distance  relation  from  Figure  31  and  the  solid  lines  represent  the 
attenuation  relations  of  the  same  slope  which  provide  best  fits  to  the 
observed  ground  motion  amplitudes  for  these  three  mine  blasts.  The  vertical 
arrows  with  annotated  values  indicate  the  average  multiplicative  factors 
separating  these  two  lines  in  each  case.  Now,  since  a  1  kt  fully  decoupled 
explosion  is  expected  to  produce  low  frequency  seismic  motions  with 
amplitudes  con^arable  to  those  expected  from  a  14  ton  (i.e.,  1000  tons/70) 
fiilly  tamped  nuclear  explosion  at  the  same  source  location,  it  might  be 
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expected  diet  the  displacement  amplitudes  from  the  1  kt  decoupled  test  would 
be  comparable  to  those  observed  from  the  S/8/92  11.9  ton  CE  test  and 
significantly  lower  than  those  observed  from  the  4/12/92(1)  27  ton  and 
4/12/92(2)  70  ton  CE  tests.  However,  it  can  be  seen  from  Figure  33  that,  in 
fact,  the  average  1  kt  nuclear  decoupled  amplitude  levels  represented  by  the 
dashed  lines  lie  well  above  the  observed  CE  amplitude  levels  in  all  three 
cases.  These  observations  are  consistent  with  the  theoretical  simulation 
results  reported  by  Barker  et  al.  (1992)  which  indicate  that  the  seismic 
coi^ling  efficiencies  of  ripple*fired,  excavation  CE  blasts  are  significantly 
lower  than  those  to  be  expected  from  contained  underground  nuclear 
explosions  of  comparable  yield.  The  comparisons  of  Figure  33  suggest  that, 
on  average,  a  CE  mine  blast  with  a  yield  of  about  75  tons  at  these  mines 
would  be  required  to  produce  the  same  low  frequency  ground  motion 
amplitudes  as  a  1  kt  fiilly  decoupled  or  14  ton  tamped  underground  nuclear 
explosion.  That  is,  the  low  frequency  seismic  coupling  efficiency  of  these 
CE  mine  blasts  is  on  the  order  of  a  factor  of  five  lower  than  that  expected 
from  a  tamped  nuclear  explosion  of  comparable  yield.  Of  course, 
significant  deviations  from  this  average  coupling  difference  can  be  expected 
in  individual  cases  due  to  variations  in  blasting  practice. 

Selected  vertical  component  ground  motions  observed  from  the  two 
4/12/92  Tymyauz  mine  blasts  are  shown  in  Figure  34  where  they  are 
coiiq>ared  with  scaled  Azgir  ground  motions  recorded  from  the  3/29/76 
cavity  decoupled  nuclear  test  in  the  same  near-regional  distance  range.  It 
can  be  seen  that  there  are  no  obvious  differences  between  these  recorded 
ground  motions  which  would  permit  an  analyst  to  visually  discriminate 
between  the  two  source  types.  Of  course,  this  conclusion  is  contingent  on 
the  unstated  assumption  that  there  are  no  major  difierences  in  regional  phase 
propagation  characteristics  between  the  Azgir  and  KIV  sites  which  might  be 
obscuring  source  related  differences  in  the  recorded  signals.  This 
assumption  will  have  to  be  tested  more  carefully  in  the  next  phase  of  this 
investigation.  In  any  case,  the  preliminary  comparisons  shown  in  Figure  34 
suggest  that  it  will  not  be  a  trivial  matter  to  discriminate  between  CE  and 
cavity  decoupled  nuclear  explosions  using  broadband,  near-regional  seismic 
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Figure  34.  Comparison  of  scaled  Azgir  1  kt  fully  decoupled  ground  motion  estimates  with  the  ground 
motions  recorded  in  the  same  near-regional  distance  range  from  CE  blasts  at  the  lymyauz  mine. 


data.  Additional  research  will  be  required  to  determine  whether  more 
detailed  and  sophisticated  comparisons  of  the  data  will  provide  a  reliable 
seismic  procedure  for  identifying  small,  cavity  decoupled  nuclear  explosions 
in  a  background  consisting  of  large  numbers  of  CE  events  of  comparable 
magnitude. 


5.  SmniiiMry  and  Conclusions 


5.1  Summary 

Successful  seismic  monitoring  of  any  eventual  Comprehensive  Test 
Ban  Treaty  will  require  the  development  of  a  reliable  capability  to  identify 
seismic  signals  from  small,  cavity  decoupled  nuclear  explosions  from  among 
the  numerous  signals  to  be  expected  from  earthquakes  and  CE  events  of 
comparable  magnitude.  The  preliminary  investigations  sunomarized  in  this 
report  have  centered  on  a  variety  of  comparative  analyses  of  observed  and 
simulated  seismic  data  corresponding  to  decoupled  nuclear  explosions  and 
different  types  of  CE  events.  In  particular,  data  recorded  from  cavity 
decoupled  nuclear  tests  conducted  in  both  the  U.S.  and  former  Soviet  Union 
have  been  used  as  points  of  reference  for  evaluating  potential  seismic 
discriminants  which  might  be  used  to  differentiate  such  sources  from  CE 
events  occurring  in  the  same  source  regions. 

Near-regional  seismic  data  recorded  from  the  U.S.  nuclear  cavity 
decoupling  experiment  STERLING  and  the  nearby  tamped  CE  test 
STERLING  HE  were  compared  in  Section  2,  where  the  time  and  frequency 
domain  characteristics  of  the  signals  were  examined  in  detail.  It  was  found 
that,  although  the  relative  spectral  compositions  of  the  P  and  S  waves 
observed  from  these  two  sources  appear  to  be  roughly  comparable,  the  S/P 
spectral  ratios  are  significantly  higher  for  STERLING  HE  than  for 
STERLING  over  a  broad  frequency  band  extending  from  at  least  1  to  10  Hz. 
Preliminary  evaluation  of  possible  causes  of  diis  difference  seems  to  rule  out 
differences  in  charge  geometry,  but  differences  in  locations  of  the  two 
sources  widiin  the  Tatum  salt  donre  structure  remains  as  a  possible 
explanation  which  needs  to  be  tested  more  rigorously  in  future  studies. 

Short-period  P  wave  data  recorded  at  NORSAR  from  Soviet  nuclear 
tests  conducted  in  salt  cavities  at  Azgir  and  from  nearby  CE  events  were 
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compared  in  Section  3  where  the  array  data  were  processed  to  obtain 
average  P  wave  spectra  for  each  selected  event.  The  results  of  this  signal 
processing  analysis  were  then  used  to  demonstrate  that  the  variation  of  the 
spectral  characteristics  of  the  recorded  P  waves  from  the  different  types  of 
CE  events  envelopes  the  spectral  characteristics  of  the  recorded  P  waves 
from  nearby  nuclear  cavity  explosions  of  comparable  magnitude.  This 
evidence  was  interpreted  as  an  indication  that  teleseismic  P  wave  data  in  the 
band  from  about  1  to  10  Hz  are  unlikely  to  provide  a  reliable  basis  for 
discriminating  between  these  two  source  types. 

Finally,  in  Section  4,  broadband  seismic  data  recorded  at  near- 
regional  Soviet  stations  from  an  Azgir  nuclear  cavity  decoupling  test  were 
theoretically  scaled  to  the  1  kt  fully  (tecoupled  level  and  con^>ared  with  data 
recorded  at  the  nearby  IRIS  station  KTV  from  well-documented  CE  events. 
Using  die  results  of  tiiese  comparisons,  it  was  shown  that  the  ground  motion 
amplitude  levels  corresponding  to  the  1  kt  fully  decoupled  nuclear  test  are, 
on  average,  a  factor  of  five  larger  than  those  observed  from  near-surface 
mine  blasts  of  conparable  energy  release,  in  agreement  with  the  results  of 
previous  finite  difrerence  simulations  of  the  relative  seismic  coupling 
effrciency  of  these  two  source  types.  Moreover,  the  broadband,  near- 
regional  ground  motion  time  history  data  recorded  from  these  events  were 
shown  to  be  highly  complex  and  variable  and  not  obviously  indicative  of 
source  type.  These  preliminary  findings  were  interpreted  as  an  indication 
that  additional,  innovative  research  will  be  required  to  define  a  reliable 
seismic  procedure  for  identifying  small  cavity  decoupled  nuclear  explosions 
using  broadband,  regional  data. 

5.2  Conclusions 

The  research  sununarized  above  supports  the  following  preliminary 
conclusions  regarding  the  seismic  identifrcation  of  small,  cavity  decoupled 
nuclear  explosions. 
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(1)  The  STERLING  HE  test  produced  seismic  signals  with  amplitude 
levels  which  are  about  twice  as  large  as  those  expected  from  a 
tamped  nuclear  explosion  of  comparable  energy  release.  It 
follows  that  a  tamped  CE  test  with  a  yield  of  about  7  tons  would 
be  expected  to  produce  ground  motion  amplitudes  equal  to  those 
expected  from  a  1  kt  fully  (tecoupled  nuclear  cavity  test. 

(2)  The  S  to  P  wave  spectral  ratios  for  STERLING  HE  arc 
significantly  higher  than  those  for  STERLING  over  a  broad 
frequency  band  extending  from  at  least  1  to  10  Hz.  This  fq)pears 
to  be  unrelated  to  differences  in  charge  geometry  between  the 
two  tests,  but  may  correlate,  ^  least  partially,  with  differences  in 
source  location  within  the  Tatum  salt  dome  structure. 

(3)  The  variations  of  the  spectral  characteristics  of  the  short*period  P 
waves  recorded  at  NORSAR  from  CE  events  near  Azgir  envelope 
the  characteristics  of  the  corresponding  P  wave  spectra  for 
nuclear  cavity  explosions  of  comparable  magnitude  at  the  Azgir 
test  site.  This  evidence  suggests  that  teleseismic  P  wave  data  in 
die  band  from  about  1  to  10  Hz  are  unlikely  to  provide  a  reliable 
basis  f<x  discriminating  between  these  two  source  types. 

(4)  Observed  seismic  data  from  nuclear  tests  can  be  theoretically 
scaled  to  other  hypothetical  source  conditions  with  reasonable 
accuracy  using  simple,  analytic  approximations  to  the  nuclear 
seismic  source  functions.  In  particular,  experience  with 
SALMON/STERLING  indicates  that  the  procedure  is  applicable 
to  the  simulation  of  ground  motions  expected  from  low  yield 
cavity  decoupled  explosions. 

(5)  Preliminary  comparisons  of  free-field  ground  motion  data 
recorded  from  a  tamped  nuclear  explosion  in  salt  at  Azgir  with 
corresponding  SALMON  free-field  data  suggest  that  the  seismic 
source  coupling  is  quite  similar  for  these  two  salt  domes. 
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(6)  Comparisons  of  the  seismic  amplitude  levels  corresponding  to  a  1 
kt  fully  decoupled  nuclear  cavity  test  at  Azgir  with  the  amplitude 
levels  observed  from  well-documented  surface  mine  blasts  at  IRIS 
station  KIV  indicate  that  the  seismic  coupling  efficiency  of  these 
ripple-fired  CE  events  is  about  a  factor  of  S  lower  than  that 
expected  fit>m  a  tamped  nuclear  explosion  of  con^arable  energy 
release.  It  follows  that  ripple-fired  mine  blasts  with  explosive 
yields  of  about  70  tons  would  be  expected  to  produce  seismic 
signals  with  amplitude  levels  equal  to  those  to  be  anticipated  from 
a  1  kt  fully  decoupled  nuclear  cavity  test. 

(7)  Comparisons  of  near-regional  ground  motion  time  histories 
corresponding  to  a  1  kt  fully  decoupled  nuclear  cavity  test  at 
Azgir  with  ground  motion  data  recorded  in  the  same  distance 
range  from  ripple-fired  mine  blasts  at  station  KIV,  reveal  no 
obvious  differences  which  would  permit  an  analyst  to  visually 
distinguish  between  the  signals  from  these  two  source  types.  It  is 
concluded  that  more  detailed  and  sophisticated  comparisons  of 
such  data  will  be  required  before  it  can  be  determined  whedier  a 
reliable  seismic  procedure  can  be  defined  for  identifying  small, 
cavity  decoupled  nuclear  explosions  in  a  background  consisting  of 
signals  from  numerous  CE  events  of  about  the  same  magnitude. 


65 


REFERENCES 


Adushkin,  V.  V.,  I.  O.  Kitov  and  D.  D.  Sultanov  (1992),  "Experimental 
Results  of  USSR  Nuclear  Explosion  Decoiq>ling  Measurements,"  Paper 
presented  at  the  14^  Annual  PL/DARPA  Seismic  Research  Syn^sium, 
PL-TR.92-2210,  ADA256711. 

Barker,  T.  G.,  K.  L.  McLaughlin  and  J.  L.  Stevens  (1992),  "Numerical 
Simulation  of  Quarry  Blast  Sources,"  S-CUBED  Technical  Report  to 
DARPA,  SSS-TR-93-13859. 

Gupta,  I.  N.,  K.  L.  McLaughlin,  R.  A.  Wagner,  T.  W.  McEliiesh,  and  M.  E. 
Marshall  (1987),  "Miscellaneous  Studies  in  Decoupling,"  Teledyne 
Geotech  Annual  Report  TGAL-87-4. 

Mueller,  R.  A.  and  J.  R.  Murphy  (1971),  "Seismic  Characteristics  of 
Underground  Nuclear  Detonations.  Part  I.  Seismic  Spectrum  Scaling," 
Bull  Seism.  Soc.  Am.,  61,  1975. 

Murphy,  J.  R.  (1969),  "Discussion  of  Pi^  by  Dr.  Springer,  M.  Denny,  J. 
Healy  and  W.  Mickey,  The  Sterling  Ejq>eriment:  Decoupling  of  Seismic 
Waves  by  a  Shot-Generated  Cavity,"  J.  Geophys.  Res.,  74,  p.  6714. 

Office  of  Technology  Assessment  (1988),  "Seismic  Verification  of  Nuclear 
Testing  Treaties,"  OTA-SC-361,  Government  Printing  Office, 
Washington,  D.C. 

Perret,  W.  R.  (196^£ij  "Free-Field  Particle  Motion  from  a  Nuclear 
Explosion  in  Sah,  Part  I,"  Sandia  Laboratory  Report  VUF-3012. 

Ferret,  W.  R.  (1986b),  "Free-Field  Ground  Motion  Study,  Project  Sterling," 
Sandia  Laboratory  Report  SC-RR-68-410. 

Rimer,  N.,  W.  Proffer,  E.  Halda  and  R.  Nilson  (1993),  "Containment 
Related  Phenomenology  From  (Themical  Kiloton,"  S-CUBED  Technical 
Report  to  DNA,  SSS-DTR-94- 14405. 

Ringdal,  F.,  G.  B.  Young  and  D.  R.  Baumgardt  (1983),  "Study  of  Detection 
and  Location  Techniques  for  Seismic  Events  Near  Azgir,  USSR:  I. 
NORSAR  Detection  and  Location  Results,"  Technical  Report  No.  1, 
SAS-TR-82-01,  ENSCO  Inc. 


66 


Riviire-Barbier,  F.  (1993),  "Status  Report  on  the  Cooperative  Research 
Program  Between  the  EME  and  the  CSS  for  Exchange  and  Analysis  of 
Seismic  Data  from  Quarry  Blasts  in  Russia  and  in  the  U.S.."  DARPA 
Center  for  Seismic  Studies  Report  C93-01. 

Springer,  D.,  M.  Deimy,  J.  Healy  and  W.  Mickey  (1968),  "The  Sterling 
Experiment:  Decoupling  of  Seismic  Waves  by  a  Shot-Generated 
Cavity,"  J.  Geophys.  Res.,  73,  pp.  5995-6011. 

Stevens,  J.  L.,  N.  Rimer,  J.  R.  Murphy,  T.  G.  Barker,  E.  Bailey,  E.  J.  Halda, 
W.  J.  Proffer,  S.  H.  Rogers  and  B.  Shkoller  (1991),  "Simulation  of 
Seismic  Signals  from  Partially  Coupled  Explosions  in  Spherical  and 
Ellipsoidal  Cavities,"  S-CUBED  Final  Report  SSS-FR-91- 12735. 


67 


Pnf.  Itaias  Ahiou 
Seiniologicil  Lab,  2S2-21 
PlviatowOeoic^lical  A  Planetaiy  Sciences 
Cafifonik  faititaiB  TeduKdogy 
P!HMktta.CA  9U2S 

PratKeiitiAki 
Cenier  for  Earth  Sciff»cfr 
Umver^  of  Soittfaem  Cafifmnia 
University  Park 
LosAng^tCA  90089-0741 

Pntf.  ^Iton  Alexander 
Department 
4(D  Dnke  BuUdmg 
The  Pennsylvania  Stale  Univeiaity 
University  Park,  PA  16802 

Pntf.  Charies  B.  Arduunbeau 
ORES 

Umveisity  of  Coloiado 
Boulder.  CO  80309 


Dr.  Thomas  C.  Bache,  Jr. 
Sdenoe  Apfdkatkms  Cocp. 

10260 Campus  Point  Drive 
San  Diego.  CA  92121  (2  copies) 


Prof.  Muawia  Barazangi 
Institnie  for  die  Study  «  die  Continent 
Cornell  Univetsity 
Ithaca.  NY  148S3 


Dr.  Jeff  Barker 

Department  of  OedogicalSdeoces 
Suae  IMvershy  <^New  York 
tt  Binghamton 
Vestal  NY  13901 

Dr.  Douglas  R.  Baumgardt 
ENSCO.Inc 
5^  Pdrt  ^lyal  Road 
Springfield.  VA  221S1-2388 


Dr.  Susan  Beck 
Dqmrtment  of  Geosciences 
Bi^din^«77 
Univetsity  of  Arizona 
T\iacon.AZ  8S721 

Dr.  T J.  Bennett 
S-CUBED 

A  Division  of  Maxwell  Ldxntories 
1 1800  Sunrise  VaUey  Drive.  Suite  1212 
Restoo.VA  22091 


Dr.  Robert  Blandford 
AFTACnr,  Center  for  Seismic  Studies 
1300  North  17th  Street 
Suite  14S0 

Arlington.  V A  22209-2308 

Dr.  Stqdien  Biatt 
ARPA/NI^O 
3701  North  Fairfax  Mve 
Arlington.  VA  22203-1714 


Dr.  Lawrence  Burdick 
K3PP.  A-025 

Scijpps  Institute  d'  Oceanogia|diy 
Univetsity  of  California,  Sw  Diego 
U  Jolla.  CA  92093 

Dr.  Robert  Bunidge 
Schlumberger-DoU  Research  Center 
CXd  Quarry  Road 
Ridgefield.  CT  06877 


Dr.  Jetty  Carter 
Cenier  for  Seismic  Studies 
1300  North  17th  Street 
Suite  1450 

Arlington.  VA  22209-2308 

Dr.  Eric  Chad 
Division  9iMl 
Sandia  Laboratory 
Albuquerque,  NM  87185 


Dr.  Martin  Qumman 
DqMtnnent  d  uedogical  Sciences 
>fitgiiiia  Pdytechnical  Institute 
21044  DerringHaU 
Blacksburg,  VA  24061 

Mr  Robert  Codcertiam 
Arms  Contid  &  Disarmament  Agency 
320  21st  Street  North  West 
Room  5741 

Washington.  DC  20451. 

Prof.  Vernon  F.  Cormier 
Dqurtment  of  Gedogy  &  Geofdiysics 
U-45,  Room  207 
.University  of  Connecticut 
Storrs,  CT  06268 

Prof.  Steven  Day 

Dqrartment  d  Gedogical  Sciences 
San  Diego  State  University 
San  Diego,  CA  92182 


1 


MwinDeniiy 
fix  Dwiaieiit  of  Baefgy 
Oflfeo  of  Anas  CcMtfiQt 
WMliiaglin.DC  20SgS 


Dr.ZotaaDer 

ENSCO.I11C. 

3400  POft  Royal  Road 
S|inagfield.VA  22151-2388 


Pndl  Adnn  Dziewonski 
fic^Eman  Lidxxatoiy.  Harvard  University 
Dept  ttfEaidi  Atmos.  &  Planetary  Sciences 
20  Oxford  Street 
Cambridge.  MA  02138 

Prof.  John  Ebel 

Department  (rf  Geology  A  Geophy^ 
BosKmCtrflene 
ChestnutW^  02167 


Eric  Fielding 

SNEEHaU 

INSTOC 

Cotndl  Univetsiar 
Ithaca.  NY  148S3 

Dr.PnrRrbas 

Institute  of  Physia  the  Earth 
Masai^  Univetsity  Brno 
Jecna^ 

612  46  Bmo.  Czech  Rqwblic 

Dr.  Mark  D.  Fisk 

Mission  RiMearch  Corpcnadon 

735  Stale  Street 

P.O.  Drawer  7 19 

Sana  Barbara,  C  A  93102 

Ftof  Hanley  Flatte 
^xdied  Sciences  Building 
Umveisity  of  California,  Santa  Cruz 
SantaCruz,CA  95064 


Prof.  Donahl  Forsyth 
DqiMutment  of  Getdogical  Sdmioes 
Brown  Umversi^ 

Providence.  RI  02912 


Dr.  ArtPraiikd 
U.S.  Geological  Survey 
922NnionalCeiiier 
Resun.VA  22092 


Dr.  Cliff  Fndich 
Institute  of  Geophysics 
8701  Nordi  Mopac 
Austin.  TX  7^59 


Dr.  Holly  Given 
IGPP,  A-025 

Scripps  Institute  oi  Oceanography 
University  of  California.  San  Diego 
U  Jolla,  C  A  92093 

Dr.  Jeffrey  W.  Given 
SAIC 

10260  Campus  Point  Drive 
San  Diego.  CA  92121 


Dr.  Dale  Qover 
Defense  Uuelligence  Agency 
ATTN:  ODT-IB 
Washington.  DC  20301 


Dan  N.  Hagedtm 
Pacific  No^west  Laboratorks 
BatldUe  Boulevard 
Richland.  WA  99352 


Dr.  James  Hannon 

Lawrence  Livermtne  National  Labtxatory 

P.O.  Box  808 

L-205 

Livermore,  CA  94550 

Prof.  David  G.  Harkrider 
Seismological  Labmtory 
Division  of  Geological  &  Planetary  Sciences 
California  Institule  of  Technology 
Pasadena.  CA  91125 

'Prof.  Danny  Harvey 
ORES 

University  of  Colorado 
Boulder,  CO  80309 


Prof.  Donald  V.  Helmberger 
Sdsmological  Laboratory 
Division  of  Geological  &  Planetary  Sciences 
California  Institule  of  Technology 
Pasadena,  CA  91125 

Prof.  Eu^ne  Herrin 

Institute  for  the  Study  of  Earth  and  Man 

Geophyacal  Laboral^ 

Southern  Methodist  University 
Dallas,  TX  75275 


2 


PiroC  Robert  B.  Hemnaiin 

Dmofineitt  of  Earth  &  Atmospheric  Sckaces 

St  Louis  Univecrity 

St  Louis.  MO  63136 


Prof.  Lane  R.  Johnson 
Seismographic  Strtkm 
Univecsity  of  California 
Berkeley.  CA  94720 


Prof.  Thomas  H.  Jordan 
Draartment  (tf  Earth.  Atmospheric  & 
Puuirtaty  Sciences 

Massachusetts  Instituie  of  Technology 
Cambridge.  MA  02139 

Prof.  Alan  Kafka 

Department  of  Geology  &  Geophysics 
Boston  Coll^ 

Chestnut  Hill.  MA  02167 


RcrtiertC.  Kemerait 
ENSCO.Inc. 

443  Pineda  Court 
Melbourne.  FL  32940 


Dr.  Kad  Koch 

Institute  for  the  Study  of  Earth  and  Man 
Get^^ical  LaboratCHy 
Soawaa  Methodist  Univetsiiy 
Danas.  Tx  73273 

Dr.  Max  Kotmtz 
U.S.  D^  of  Enetgy/DP  3 
Fonestal  Building 
1000  Indqiendaice  Avenue 
Wariiington,DC  20383 

Dr.  Richard  LaCoss 

MIT  Lincoln  Laboratory.  M*200B 

P.O.  Box  73 

Lexington.  MA  02173-0073 


Dr.  Fred  K.  Lamb 

Universtty  of  Illinois  at  Urbana-Champaign 
Dnartment  of  Phy^ 

1 1 10  West  Green  Street 
Urbana.lL  61801 

Pnrf.  Quries  A.  Langsttm 
Geosdences  Dnartment 
403DeikBBuil^ 

The  Pranqdvania  State  Univnsity 
University  Parit  PA  16802  3 


Jim  Lawson.  Chief  Geophysicist 
Oklahoma  Geological  Survey 
Oklahoma  Geophysical  Observatory 
P.O.  Box  8 

Leonard.  OK  74043-0008 

Prof.  Thorne  Lay 

Institute  of  Tectonics 

Earth  Science  Board 

University  of  California,  Santa  Cruz 

SanUCruz.CA  93064 

Dr.  William  Leith 
U.S.  Geological  Survey 
Mail  Stop  9^8 
Reston,  VA  22092 


Mr.  James  F.  Lewkowicz 
'Phillips  Laboratory/GPffl 
29  Ramkdph  Road 

Hanscom  AFB,  MA  01731-3010(  2  copies) 


Mr.  Alfied  Lieberman 

ACDA/VI-OA  State  DqMilment  Building 

Ro(Hn3726 

320-21st  Street.  NW 

Washington.  DC  20451 

Prof.  L.  Timothy  Long 
School  of  Geophysical  Sciences 
Georgia  Institute  of  Technology 
Atlanta.  GA  30332 


Dr.  Randol|di  Martin,  III 
New  England  Research,  Inc. 

76  CMcott  Drive 

White  River  Junction,  VT  05001 


Dr.  Robert  Masse 
Denver  Federal  Building 
Box  25046,  Mail  Stop  961 
Denvn-.CO  80225 


Dr.  Gary  McCartor 
Dq>aiti^t  of  niysics 
Southern  Mohodist  Univoaty 
Dallas.  TX  75275 


Prof.  Thomas  V.  McEvilly 
Seismographic  Statioi 
University  of  California 
Berkeley.  CA  94720 


Gt.  AftMcCter 
U.S.  Geological  Survey 
Midi  Stop  9^ 

U.S.  Oedogical  Snrv^ 

Menlo  Park.  CA  9402S 

Dr.  Kekh  L.  McLaughlin 
S<XJBED 

A  Diviaton  of  Maxwell  LidKxiUory 
P.O.  Box  1620 
U  Jolla.  CA  92038-1620 

Steuben  Miller  &  Dr.  Alexandn*  Florence 

SRikiinnational 

333  Ravenswood  Avenue 

Box  AF 116 

MenloPaik.CA  94025-3493 

Prof.  Bernard  Minster 
IGPP.  A-025 

Scripps  Institute  of  Ooeanogr^y 
Univeisity  of  CaBfomia.  Sim  Diego 
Uldlla.CA  92093 

Prof.  Brian  J.  Mitchell 

Department  of  Bar^  &  Atmospheric  Sciences 

St.  Louis  University 

Sl  Louis.  MO  63156 


Mr.  Jade  Murphy 
S-CUBED 

A  Divuioo  d’Maxwdl  Laboratory 
11800  Sunrise  Valley  Drive.  Suite  1212 
Reston.  VA  22091  (2,  Cojpks) 

Dr.  Keith  K.  Nakanishi 

Lawrence  Livermore  National  Laboratory 

L-Q25 

P.O.  Box  808 
Livermore.  CA  94550 

Prof.  Jdm  A.  Orcutt 
IGPP.  A-025 

Scripps  Institute  of  Oceanography 
University  of  Oalifornia.  San  Diego 
U  Jolla.  C  A  92093 

Prof.  Jeffiey  Park 
Kline  Gedogy  Labenatory 
P.O.  Box  6^ 

New  Haven.  CT  06511-8130 


Dr.  Howard  Pamm 

Lawrence  Livermore  National  Laboratory 
L4125 

P.O.  Box  808 
Livermore.  CA  94550 


Dr.  Frank  Pilotte 
HQAFTACmr 
.  IC^  SouA  Highway  AlA 
Patrick  AFB.FL  32925-3002 


Dr.  Jay  J.  Pulti 
Radix  Systems,  Inc. 

201  Perry  Parkway 
Gaithersburg,  MD  20877 


Dr.  Robert  Remlce 
ATTN:  PCTVTD 
Field  Conunand 
Defense  Nudear  Agency 
KirtlandAFB.NM  87115 

Prof.  Paul  G.  Richards 
Lamont-Doherty  Geological  Observatory 
oi  Columbia  University 
Palisades.  NY  10964 


.Mr.  Wilmer  Rivets 
'  Teledyne  Geotech 
314  Montgomery  Street 
Alexandria.  VA  22314 


Dr.  Alan  S.  RyaU.  Jr. 
ARPA/NMRO 
3701  North  Fairfax  Drive 
Arlington.  VA  22203-1714 


Dr.  Ridiatd  Sailor 
TASC,  Inc. 

55  Walkers  Brook  Drive 
Reading,  MA  01867 


Prof.  Charles  G.  Sammis 
Center  for  Earth  Sciences 
University  of  Southern  Califrmiia 
University  Park 
Los  Angeles.  CA  90089-0741 

Prof.  Christopher  H.  Scholz 
Lamont-Dobetty  Geological  Observatory 
of  Cdumbia  University 
Palisades.  NY  10964 


Dr.  Susan  Schwartz 
Institute  of  Tectonics 
1 156  High  Street 
Santa  Cruz,  CA  95064 


4 


SeoMMy  of  die  Air  Force 
(SAFRD) 

Wa8iiiiiiloii.DC  20330 


Office  of  die  Secmaiy  df  D^Bose 
DDRAE 

Wa8iiiiiitaii.DC  20330 


Thomas  J.  Sereno.  Jr. 

Scknoe  ^iplicadon  Inti  C(»p. 
10260  Campus  Point  Drive 
SanDieio.CA  921?! 


Dr.  Nfidmel  Shorn 
D^iease  Nuclear  Agency^PSS 
6801  Teleiiaph  Road 
Alexandria.  VA  22310 


Dr.  Robert  Shumway 
Umveisity  of  California  Davis 
Divisimi  of  Statistics 
Davis.  CA  9S616 


Dr.  Matthew  Sibol 
VhiiniaTech 
Sdimok>ilca]  Observatory 
4044  Derring  Hall 
Blacksburg.  VA  24061-0420 

Prc^.  David  O.  Simpstm 
BUS.  Inc. 

1616  Noth  Fort  Myer  Drive 
Suite  IQSO 
Arlington.  VA  22209 

Donald  L.  Sfninger 

Lawrence  Livermore  Nadcmal  Laboratory 
L4I2S 

P.O.  Box  808 
Livermore.  C  A  94SS0 

Dr.  Jeffrey  Stevens 
S-CUBED 

ADtvision  (rf* Maxwell  Laboratory 
P.O.  Box  1620 
La  Jolla.  CA  92038-1620 

Lt  Col.  Bm  Stobie 
ATTN:  AFOSR/NL 
110  Duncan  Avemie 
Bolling  AFB 

Washington.  DC  30332-0001 


Brian  Stuap 

Los  Alauos  National  Laboratory 
EES-3.  Mall  Stop  C335 
Lob  Alaaos.  NM  87545 


Prof.  Jmemiah  Sullivan 
University  of  Illinois  at  Urbana-Champaign 
•Departmmt  of  Phyrics 
1 1  to  West  Green  Street 
Urbana,  IL  61801 

Prof.  L.  Sykes 

Lamont-Doherty  Geological  Observatory 
of  Columbia  University 
Palisades.  NY  10964 


Dr.  David  Taylor 
ENSCO,  Inc. 

445  Pine^  Court 
Melbourne.  FL  32940 


Dr.  Steven  R.  Taylor 

Los  Alamos  National  Laboratory 

P.O.  Box  1663 

Mail  Stop  C335 

Los  Alamos.  NM  87545 

.Prof.  Clifford  Thurber 
University  ol  V^sconsin-Madison 
Department  of  Geology  &  Geqihysics 
1215  West  Dayton  Strett 
Madison.  WS  53706 

Prof.  M.Nafi  Toksoz 
Earth  Resources  Lab 
Massachusetts  Institute  of  Technology 
42  Carleton  Street 
Cambridge,  MA  02142 

Dr.  Larry  Turnbull 
OA-OSWR/NED 
Washington,  DC  20505 


Dr.  Gregory  van  der  Vink 
IRIS,  Inc. 

1616  Neath  Fort  Myo*  Drive 
Suite  1050 

Ariington,VA  22209 

Dr.  Karl  Veith 
EG&G 

5211  Audi  Road 
Suite  240 

Suitland.MD  20746 


5 


PratTenyC  Wallace 

Umveta^  (tf  Aiizana 
‘niscoo.AZ  8S721 

Dr.  Thomas  Weaver 

Los  Alainos  Natioiial  Laboratory 

P.O.  Box  1663 

Mail  Stop  C335 

Los  Alamos,  NM  87S4S 

Dr.  William  Wortman 
Mtsskm  Reseavdi  Corponttkm 
8360  Ondeibed  Road 
SoileTOO 

Newington,  VA  22122 

Prof.  Prancis  T.  Wu 
Department  of  Gecriogical  Sciences 
Stale  (laiversity  of  York 
atBin^hamtoa 
VestaUfY  13901 

Prc^Ru-Shan  Wtt 

Untveisiiy  of  California,  Santa  Cruz 
Earth  Sdoioes  Department 
Santa  Cniz 
,CA  93064 

ARPA,OASB/Ubniy 
3701  NbrA  FidtfEa  IXive 
Arlington,  VA  22203-1714 


HQI»4A 

ATTN:  Techmcal  library 
Washington.  DC  20305 


Ddoise  Intelligenoe  Agency 

Ditedonae  for  Scientific  &  Technical  hitdligenoe 

ATTN:  DTIB 

Washington.  DC  20340-6138 


Defense  Technical  Infbrmatkm  Center 
Camenm  Station 

Alexandria,  VA  22314  (2  Copies) 


TACIEC 

Baaidle  Memorial  histitute 
305  King  Avemie 

Odombra,  CXI  43201  (Final  Rqxrit) 


Phillms  Laboratory 

ATTN:  XPG 

29  Randolph  R(»d 

Hanscom  AFB.  MA  01731-3010 


Phillips  Laboratory 

ATTN:  GPE 

29  Randolph  Road 

Hanscom  AFB.  MA  01731-3010 


Phillips  Laboratory 
ATTN:  TSML 
-5  Wrij^t  Street 

Hanscom  AFB.  MA  01731-3004 


Phillips  Laboratory 
ATTN:  PL/SUL 
3530  Aberdeen  Ave  SE 
Kirtland.NM  871 17-5776  (2  copies) 


Dr.  Michel  Bouchon 
I.R.I.G.M.-B.P.  68 
38402  Sl  Martin  DUeres 
CCdex.  FRANCE 


Dr.  Michel  Campilk) 
Observatoire  de  Grerwble 
I.R.I.G.M.-B.P.  53 
38041  (}renoble.  FRANCE 


.Dr.  Kin  Yip  Chun 
Geophysics  Division 
Physics  Dqrartment 
University  of  Toronto 
Ontario,  CANADA 

Prof.  Hans-Peter  Har^ 
Institute  for  Cieoidiymc 
Ruhr  University^ochum 
P.O.  Box  102148 
4630  Bochum  1.  GERMANY 

Prof.  Eystein  Husebye 
NTNF/NORSAR 
P.O.  Box  51 

N-2007  Kjeller,  NORWAY 


David  Jepsen 

Acting  Head,  Nuclear  Mcmitoiing  Sectitxi 
Bureau  of  Mineral  Resources 
Geology  and  Geophysics 
G.P.O.Box378,  Canberra.  AUSTRALIA 


Ms.  Eva  Johaniusson 
Senior  Researdi  Officer 
PQA 

S-172  90  Sundbyberg,  SWEDEN 


*  Dr.FblerManliaU 
I^nocurement  Executive 
Ministry  ofDefense 
^  Bladmest,  Biimpton 

Reading  FG7-FRS.  UNITED  KINGDOM 

Dr.  Bernard  Massinon.  Dr.  Pierre  Mechler 
Societe  Radiomana 
27  rue  Claude  Bernard 
75005  Paris.  FRANCE  (2  Copies) 


Dr.  Svein  Mykkeltvdt 
NTNT/NORSAR 
P.O.  Box  51 

N-2007  Kjeller.  NORWAY  (3  Copies) 


Prof.  Keith  Priestley 
Univeraty  of  CamMdge 
Bullard  Ud»,  Dept  of  Earth  Sciences 
Madingl^  Rise,  Madingley  Road 
Cambndge  CB3  OEZ.  ENGLAND 


Dr.  Jorg  Schlinenhardt 

Federal  Institute  for  Geosciences  &  Natl  Res. 
Postfach  510153 

D-30631  Hannover ,  GERMANY 


Dr.  Johaimes  Schweitzer 
Institute  of  Geophysics 
Ruhr  University/Bochum 
P.O.  Box  1102148 
4360  Bochum  1.  GERMANY 

Trust  &Vmfy 

VERHC 

Carrara  House 

20  Embankment  Place 

London  WC2N  6NN.  ENGLAND 


7 


